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Abstract 

We study properties of the compact 4D U(l) lattice gauge theory with monopoles 
removed. Employing Monte Carlo simulations we calculate correlators of scalar, 
vector and tensor operators at zero and nonzero momenta p . We confirm that 
the theory without monopoles has no phase transition, at least, in the interval 
< (3 < 2 . There the photon becomes massless and fits the lattice free field theory 
dispersion relation very well. The energies of the ++ , l + ~ and 2 ++ states show 
a rather weak dependence on the coupling in the interval of (3 investigated, and 
their ratios are practically constant. We show also a further modification of the 
theory suppressing the negative plaquettes to improve drastically the overlap with 
the lowest states (at least, for J = 1). 
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1 Introduction 



Despite the big efforts invested last years into the nonperturbative study of QED 
(for reviews see, for instance, [I], |2|) we still have no clear understanding of this 
theory. For the standard compact formulation there is a strong indication for a 
first order phase transition, thus leaving no room for the continuum limit. The 
noncompact formulation seemed to be more promising because the (chiral) phase 
transition turned out to be of second order. However, the noncompactness is realized 
only in the pure gauge part of the action, while in the fermionic part a compact 
U(l) phase remains. This part of the action will eventually induce the compact 
term back into the action, thus living us with the question of selfconsistency of 
the (so-called) noncompact theory. Moreover, the compact part of the action is 
responsible for the appearance of monopoles, which are lattice artifacts and are 
supposed to be at least in part the driving force of the chiral phase transition. If 
the latter is correct then the chiral phase transition in the noncompact QED might 
not have a continuum analogue ||. 

The problem of lattice artifacts should be taken very seriously in lattice cal- 
culations. Within the Monte Carlo approach, in order to calculate the average of 
any field operator one generates a sequence of (equilibrium) gauge field configura- 
tions: {Uil^ k }, {Uj?lk}, ■ ■ ■ with the statistical weight P(Ui) ~ exp [— S(Ui)] , 
where S(Ui) is the Euclidean lattice discretized action, and averages over all 
of them. The crucial question in this approach is, whether in a certain range of 
coupling(s) these field configurations can provide an adequate representation of con- 
tinuum physics. It is a well-known fact that artifacts can change the structure of 
a theory drastically. One example of the influence of lattice artifacts is provided 
by the study of topological properties of gauge theories. Small-scale fluctuations, 
i.e., the fluctuations living on the scale size of, say, one lattice spacing, and carry- 
ing a non-trivial topological charge (dislocations) can lead to the divergence of the 
topological susceptibility in the continuum limit (at least for the geometric defini- 
tion of the topological charge). A special source of troubles are field configurations 
containing negative plaquettes: |Tr£/p ~ —1 (for SU(2) ) J|, |, ||. A similar 
situation takes place for the CP n_1 model with n < 3 0. The fact that sin- 
gular configurations may be important in quantum field theory was demonstrated 
also for the three-dimensional non-linear 0(3) model || and for the compact 
Sine-Gordon model || . It was speculated that point-like singularities (topological 
defects) are the sources of a topological anomaly, and a new coupling constant con- 
trols their behaviour. This 'hidden' coupling constant without making important 
contributions to the short-distance behaviour becomes crucial in understanding the 
infrared structure (long-distance behaviour) of the theory. Rough fluctuations can 



also exhibit themselves in the form of Zn— strings and Zn— monopoles ||iO| - [l5|1 



( Zm~ artifacts). In the paper |L6| the influence of small-size ^-artifacts (nega- 
tive plaquettes) on 'physical' values was studied for pure gauge SU{2) theory. It 
was shown that these Z^— artifacts can strongly influence large-scale objects, e.g., 
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Wilson loops and their ratios. In this connection it could be instructive to note 
the following. The number of artifacts (e.g., negative plaquettes) N_ decreases 
exponentially with increasing (3 . But this in itself does not mean that they be- 
come unimportant at large (3 . They can still strongly influence any lattice average 
(O(U))(0) (say, Wilson loops) decaying faster then iV_(/3) with increasing (3 , 
and, therefore, even a small admixture of these rapidly varying fields can become 
competitive. 

The question of contamination from lattice artifacts should be addressed also 
to the pure gauge U(l) lattice theory which is expected to be the free photon's 
theory in the continuum limit. During last years this theory was rather well studied 
(for a noncomplete list of references see, for instance, |!7| - p8i ). On the lattice the 
pure gauge U(l) theory in four dimensions with Wilson action has two phases 
separated by a first order phase transition. In the strong coupling phase monopoles 
(i.e., artifacts) are supposed to form the condensate which causes the confinement 
of electric charges. At the phase transition point (3 C this condensate dissolves, and 
the transition to the Coulomb phase occurs. This theory has a nontrivial spectrum 
which differs within these phases. Its ++ , l + ~ and 2 ++ states were studied in 



To recover the continuum limit of this theory it is important to study the role of 
the underlying singular configurations to be separated from the 'real' physics. One 
can modify the lattice action in such a way that the unphysical short-range fluctu- 
ations (monopoles, . . .) become suppressed leaving the 'physical' fields untouched. 
The question of interest is then whether this modified theory can serve as a better 
lattice approximation to that of noninteracting photons in the continuum. To study 
this question in detail is the goal of the present work. 

The paper is organized as follows. First we discuss the lattice monopoles and the 
modification of the action. Section 3 is devoted to the discussion of correlators and 
'glueball' wave functions. Section 4 presents the results of Monte Carlo calculations 
and their analysis. The last section contains a discussion of the results as well as 
some speculations. 

2 Monopoles and modified lattice action 

The standard Wilson action (WA) for the pure gauge U(l) theory is 



where f3 = l/gl are = l/g 2 , and U x = U Xfl = exp(i9 x . „), 9 X -^ G (— 7r, 7r] are the 
field variables defined on the links I = (x, fj,) . Plaquette angles 6 P = 9 X .^ U are 
given by 



E0I. 



S w (Ui) = (3 ■ ^2 (1 - cos6 x . ifU/ ) , 



(1) 



e 



&x: a + 0, 



x+fi; v 



-6, 



x-\-i>\ fi 



-e. 



X] V 



(2) 
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To extract magnetic monopoles the plaquette angle 9p is split up 

9p = 9 P + 27m P , -n<9p<7r, n P = 0,±l,±2, (3) 

where 9p describes the electromagnetic flux through the plaquette and rip is the 
number of Dirac strings passing through it. The net number of Dirac strings going 
out of an elementary 3D cube determines the monopole charge within this cube as 
follows. The monopole current out of the cube c n<n , labeled by the dual link (n, /x), 
is defined to be |1| 

K n u — £uvp<T A u @n,pa j (4) 
47T 

where the lattice derivative A u is defined by A u f n = f n+ c — f n . Then we have 

2nK n , a = Y,0p = Y.( d P- 2 ™p) = -27T^n P , (5) 
p p p 

where the sums are over oriented plaquettes P in the surface of the cube c„ M . 

£A^ = 0, (6) 

which means that the monopole currents form closed loops on the dual lattice. 
To suppress monopoles one can use the modified action (MA) 

Sma = Sw + A k ■ 1 K c 1 , (7) 

c 

where the parameter \k{(3) plays the role of a chemical potential for monopoles. 
In our calculations we have chosen A^ = oo . The partition function is 

Zma = J U dU i Il^ciO • exp[ - S w (U t ) } (8) 

links c 

with periodic boundary conditions, and the average of any field functional 0(U) 
is defined as 

(0(U)) MA = Z M \ ■ J XldUi 0{U) ■ n^ c; o • exp[ - S w (U t ) ] . (9) 

links c 



The action Sma was used in p2| , |23| in order to investigate the behaviour of the 
average plaquette (P) . 

For some (3- values we employed a further modification of the action (MAI) 
suppressing all lattice artifacts characterized by negative plaquette values 

Smai = Sma + A P ■ ^(l - sign(cos6» P )) , A P = oo . (10) 
p 

Note that both modifications (MA and MAI) do not influence the formal continuum 
limit and do not change any perturbative aspects. 



4 



3 Correlators and boson wave functions 



Energies of the excited states ('glueballs') can be calculated by measuring the 
asymptotic behaviour of the correlation functions of operators $ with the proper 
J PC quantum numbers and momenta p . 

The corresponding correlation functions T(r) on a L t ■ L z s lattice are defined 
as follows 



r(r) = ($*(t + r)-$(t)) c 

i n i Lt — 1 

= — ,ErE K*(*© T ) • Mt) - e t) • ¥(*)] , (n) 

n ~~ 1 i=l ^ t=0 



where n is the number of measurements, t © r = (i + r) mod(L() and 

_ 1 n 

m = -EHt) ■ (12) 

n r: 

This definition provides us with unbiased estimators for T(r) . 

The correlators defined above are the superpositions of the exponents corre- 
sponding to different energy levels, e.g., 



T(t) 



A ■ 



e -E-r + e -E-(N t -r) 



+ ..., (13) 



r(o) 

where dots correspond to higher state contributions and 

A = fL-\( jPC \* jP ° °)l 2 <1- ( 14 ) 

In principle the lowest state energy can be defined from the ratio of correlators 
at large enough r (but r C L f /2 ) : 

T(t) 

£ e // = -m f ^ T y, (15) 

and E = E e ff if E-{r — 1) ^> 1 . In practice only for comparatively small 
values of energy (E ~ 1 — 2) we can use the separation distance up to r = 5 . 
For larger values of energy the correlators decay very quickly, the signal-to-noise 
ratio becomes too bad, and at E^A it is possible to use only r = and r = 1 
values with reasonable statistical errors. In this case the effective energy E e ff(r) 
in eq.(|T5D is nothing but an estimate of the true energy E . The question how 
close E e ff(r) is to the true value E is the question of contamination from higher 
energy states. The larger the overlap of our 'glueball' wave function with the lowest 
lying 'glueball' state, the closer A is to 1 , the less the contamination from higher 
energy states. The common way to suppress the higher states contribution is to 
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consider 'smeared' wave functions which are linear combinations of large numbers 
of Wilson loops and have a large overlap with low lying states f3fj, |3l], 32 1. 



In our work we used the simplest plaquette wave functions which are known to 
have comparatively small overlaps with the lowest glueball states (at least for the 
standard Wilson action). For us this is not the point of great importance because 
our main goal is the comparison of predictions for the different actions. Also the 
plaquette wave functions are proved to work very well in the case of comparatively 
small energies (where we can use the larger values of r). We shall come back to 
this question lateron, when discussing the influence of ^-artifacts. 

To construct the one-plaquette operators we followed the standard recipe de- 
scribed elsewhere (for a review see, for instance, f29|). The first step is to construct 
local operators $(i;r) whose p = 0— parts have the desired J PC properties. 
Real parts of Wilson loop operators have C— parity C — +1 and imaginary parts 
have C— parity C — — 1. For the ++ state we first sum up the real parts of all 
six plaquettes belonging to the smallest space-like cube 

Ht;r) = \ £cos(0 P J , (16) 
k=l 

where f is the centre of the cube. This construction is invariant under 7r/2 and 
7i rotations about all three axes in the 3D space passing through its centre. It is 
easy to show that the zero momentum operator 



$(t)=X>(t;f) (17) 

r 

corresponds to J PC = ++ , 4 ++ , . . . exitations. We assume that the state with 
J = has the minimal mass, so, we shall refer to it as the ++ state. 

The Fourier transform with momentum p gives us non-zero momentum oper- 
ators: 



5>-**-*(*;t0 (18) 

r 

For constructing a tensor operator one has to take plaquettes parallel to the 
quantization axis (say, OY) and subtract from it the same but rotated by tt/2 
about OY . After summing through the three-dimensional space we obtain the 
operator corresponding to J PC = 2 ++ , 4 ++ , . . . excitations. As far as it is 
commonly assumed that the state with J = 2 has the minimal mass it is refered 
as the 2 ++ state with spin projection ± 2 onto the quantization axis. In this 
case the momentum p in eq.(|I^) must be chosen parallel to the quantization axis 
to avoid the admixture of a scalar state. 

Finally, summing up the imaginary parts of plaquette operators parallel to the 
quantization axis results in 3 + ~, ... states (or just l +_ state). Comparing 

with the photon it has the wrong parity at zero momentum. But at p ^ it has a 
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non-zero overlap with the 1 (photon) state [20 . This fact can be easily checked 
in the naive continuum limit. 



4 Monte Carlo calculations 

Our calculations were made on a L t ■ = 12 • 6 3 lattice at 0.1 < (3 < 2 . For each 
(3 value our statistics reaches from 40.000 to 50.000 sweeps, and 2000 sweeps 
without measurements were done to reach equilibrium. For calculating the mean 



square errors and biases of our estimators we used the jackknife method |j5[ [3(j . 

We calculated autocorrelation functions C auto \n sweep ) for all our observables 
(plaquettes P Sji , correlators T(t;p) , etc.) and defined the decorrelation length 
ridec. as a minimal number of sweeps such that 

C auto in dec .)<l-C auto i0) . (19) 
5 

For the standard Wilson action ridec. gets a sharp peak at (3 = (3 C ~ 1 with 
^cZec.~ 100 , while for the modified actions this value is comparatively small [rtdec. ~ 
4 — 6) and has a very weak dependence on f3 in the chosen interval. 

In Fig.l we show the dependence of the average plaquette (P) on (3 for WA 
and MA. For the modified action this dependence is very smooth and shows no 



trace of crossover or phase transition in agreement with the calculations of [22, 23]. 

Another important parameter which we care about is the average number of 
negative plaquettes (AL) , i.e., plaquettes P x;fJjU with 



cos(^) < . (20) 

In Fig. 2 we show the dependence of its density p neg . = (N_)/N P on (3 for WA and 
MA. It is interesting to note that the suppression of monopoles decreases the density 
of negative plaquettes at (3^1.1 , while at larger (3 values p neg , is approximately 
the same for both the actions. 

As an example in Fig. 3 we show the r dependence of the vector state correlator 
T(r) at nonzero momentum p = (0,1,0) ■ 7r/3a at /3 — 0.3 and (3 = 1.0 
for WA and MA. For the modified action the correlator does not depend on the 
coupling giving an indication for the decoupling of the photon. Because the energy 
is comparatively small these correlators have reasonable errors up to r = 5 . On 
the contrary, for the Wilson action at (3 — 0.3 we recover the expected strong 
fall-off of the correlator. 

In order to extract masses from the effective energy E e ff(p;r) we need a 
dispersion relation (DR). In the free field boson case the propagator on the lattice 
provides us with the DR 
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2 2^2' v ' 

which at small enough spacing and momenta gives us the continuum DR: 

E 2 = m 2 + p 2 . (22) 

In general the connection between the momentum and the energy can be very 
nontrivial on lattice, and a DR which we can expect is of the form (see, e.g., [ |33|j ) 



sinh 2 ^ = zMJ PC )-^h 2 ^ + z 2 (frJ PC ).j:sin 2 ^ 

+zM J pc ) ■ £ - 2 ■ sin 2 °f ■ sin 2 °f + . . . . (23) 
We believe that in the continuum limit 



Zi{(3]J PC ) -> 1 at z = l,2; 

J PC ) -> at i > 3 . (24) 



The second line in eq . (p4|) provides the necessary condition for the restoration of 
rotational invariance (at small enough p ), and the first one yields the necessary 
condition for the restoration of Lorentz invariance. 

In Figs. 4 we show the effective energy of the vector boson E e ff(p;r) as a 
function of (3 for Wilson action (Fig. 4a) and for the modified action (Fig.4b). 
Different symbols correspond to different momenta p as indicated in the figures. 
Broken lines correspond to the free field theory dispersion relation eq. fl2~T|) with zero 
mass and different momenta p . In the theory with the MA a massless photon exists 
in the whole f3 region we have chosen, and the free field theory dispersion relation 
(E21p works astonishingly well. Therefore, it provides a very convincing evidence 
that in the theory without monopoles a massless photon exists even at very small 
f3 . In the theory with the WA a massless photon exists in the Coulomb phase as 
already observed in pDI (see also [28]) while below the critical point the photon 



acquires a nonzero mass. 

Another interesting observation is that the DR (|21 ) works very well for both 
the actions at even for r = 1 . This fact shows that the overlap coefficient 

A(/3) is close enough to unity in this phase, and the contamination from higher 
states becomes negligible with increasing j3 . 

Analogously to Figs. 4 in Figs. 5 we show the effective energy E e ff(p; r = 1) of 
the ++ boson as a function of (3 for the Wilson action (Fig. 5a) and modified 
action (Fig. 5b). For the MA the dependence on the coupling is very smooth, and at 
increasing (3 the energies for both actions tend to be rather close to each other. It 
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is difficult to check the validity of the free field theory DR with our data because the 
energies are comparatively big, and their dependence on the momentum is weak. 

Similar observations can be made for the tensor states. In Figs. 6 we show the 
effective energy E e ff(p;r = 1) of the 2 ++ boson as a function of j3 for the 
Wilson action (Fig.6a) and the modified action (Fig. 6b). 

Masses m 7 we calculated are in units of the (unknown) spacing a{(3) : 

= m], h _- a(/3) ■ 7 = ++ ,... ; (25) 

but their ratio 

m 2++ 1 m 0++ = mf h + / (26) 

must not depend on the coupling in the continuum limit. Comparing the data 
for scalar and tensor states it is easy to see that this ratio has a very weak (3 
dependence in the whole interval and is rather close to unity. 

As it was shown above the simple one-plaquette wave functions work very well 
for r > 2 (at least, for the vector state). This means that the overlap coefficient 
A(P) which can be defined as 

A ^ = w^i ■ (27) 

is not too far from unity. Fig. 7a shows its dependence on j3 for the vector state at 
p = (0; 1; 0) • n/3a . For MA (open circles) A ~ 0.65 at small (3 , and at /3 ~ 1 it 
becomes close to 1 . It is tempting to assume that there is some correlation between 
this dependence and the (3— dependence of the number of negative plaquettes (com- 
pare with Fig. 2). To check this assumption we used another modified action (MAI) 
with the full suppression of negative plaquettes (|TUp. Calculations were made at 
j3 = 0.3 and (3 = 0.5 , and the corresponding values of A{(3) are shown in Fig. 7a 
with full circles. In Fig. 7b we show the effective energy E^p^p; r = 1) at (3 — 0.3 
as a function of p 2 = J2i=i 4 sin 2 (api/2) . Open circles correspond to MA, and 
full circles correspond to MAI. The solid line corresponds to the dispersion rela- 
tion (pTF). Therefore, after suppressing negative plaquettes the agreement between 



the r = 1 effective energy and the dispersion relation (|21| ) becomes practically 
perfect. So, we conclude that the Z 2 -artifacts are, at least partly, responsible for 
not good enough overlap with the higher states. It is worthwhile to note here that 
the suppression of negative plaquettes does not influence the energy of the photon 
while the energies of the massive scalar, vector and tensor states tend to increase. 



5 Conclusions 

We have studied the properties of the compact four-dimensional U{1) lattice gauge 
theory without monopoles. Employing Monte Carlo simulations we calculated cor- 
relators of scalar, vector and tensor operators at zero and nonzero momenta p on 
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a 12 • 6 3 lattice at 0.1 < f3 < 2 . We confirm that the theory without monopoles 
has no phase transition, at least, in the interval of f3 studied. In the modified 
theory a massless photon exists in the whole interval of f3 , fitting the lattice free 
field theory dispersion relation very well. The energies of ++ , l + ~ and 2 ++ 
states show a rather weak dependence on the coupling in the chosen interval of (3 , 
and their ratios are practically constant. We show also that other lattice artifacts 
- negative plaquettes - are, at least partly, responsible for the overlap in the strong 
coupling region. The further modification of the theory by suppressing the nega- 
tive plaquettes drastically improves the overlap with the lowest states (at least, for 

We believe that the problem of artifacts in lattice calculations deserves much 
more detailed and thorough study, and still our understanding of pure gauge U(l) 
theory is far from being complete. Even after suppressing monopoles (and negative 
plaquettes) there is still some residual interaction which occurs in the appearence 
of massive states (0 ++ , . . .). One can speculate that there are other kinds of lattice 
artifacts which are responsible for this interaction, and Dirac strings are among the 
possible candidates for a further study. 

We conclude that the modified theory without lattice artifacts (monopoles, neg- 
ative plaquettes, etc.) should be a more reliable candidate for the QED calculations 
with fermions in the strong coupling area. 
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Figure captions 



Fig.l The plaquette value (P) for WA (squares) and MA (circles) as a function 
of p . 

Fig. 2 The density of negative plaquettes p neg , for WA (squares) and MA (circles) 
vs. p . 

Fig. 3 Correlator T(r;p) for the vector state at P — 0.3 and P = 1.0 for WA 
and MA. p = k- 2ir/aL s 

Fig. 4 The effective energy of the vector state E e ff(P; r) at different p = k-2n /aL s 
as a function of P for WA (a) and for MA (b). Different symbols correspond 
to different momenta p shown in the picture. Broken lines correspond to the 
dispersion relation ( pT|) with zero mass. 

Fig. 5 The effective energies E e ff(p; r = 1) for the scalar state as a function of ft 
for WA (a) and for MA (b). 

Fig. 6 The effective energies E e ff(p; r = 1) of the tensor state as a function of P 
for WA (a) and for MA (b). 

Fig.7 A(P) = T 2 (r = l)/r(r = 2) vs. p for the vector state for MA and MAI 
(a); The effective energy E e ff(p;r = 1) of the vector state at P = 0.3 as a 
function of p 2 = J2i=i 4 sin 2 (api/2) for MA and MAI. The momenta are the 
same as in Fig. 4. The solid line corresponds to the dispersion relation ([H]) (b). 
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